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Abstract 

In this paper, we find all integers x such that x^ — 1 has only 
prime factors smaller than 100. This gives some interesting numerical 
corollaries. For example, for any positive integer n we can find the 

^D ' largest positive integer x such that all prime factors of each of x, a; + 

^N , 1, . . . , a; + n are less than 100. 

>>■ 
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■ 1 Introduction 

H, 

^ ' For any integer n we let P{n) be the largest prime factor of n with the 

r~| . convention i-*(0) = P(ibl) = 1. Our main result in this paper is the deter- 

j^ \ ruination of all integers x satisfying the inequality 

P(a;2 - 1) < 100. (1) 

We also give some interesting corollaries to our main result. For simplicity, 
we will consider only solutions with positive values of x. 
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Before stating the results, some history. In 1964, Lehmer [B] found all 
positive integer solutions x to the inequality P{x{x + 1)) < 41. Notice 
\f^ . that this amounts to finding all odd positive integers y = 2x + 1 such that 

^D \ P{y'^ — 1) < 41. There are 869 such solutions. In ^, the first author found 

^~^ ■ all positive integer solutions of the similar looking equation P(x^ + 1) < 100. 

There are 156 of them. 

In both [B] and [7], the method of attack on this question was the follow- 
ing. Assume that j; is a positive integers such that P{x'^ ±1) < K, where 
p^ ', K = Al if the sign is — and K = 100 if the sign is +. Then we can write 

x^±l = dy\ (2) 

where d is squarefree, and P{dy) < K. When the sign is — , then there are 
13 possible primes p < 41 which can participate in the factorization of d. 
When the sign is +, then p < 100 but either p = 2 or p is a prime which is 
congruent to 1 modulo 4. There is a total of 12 such primes. In both cases, 
we can write equation ([2]) in the form 

x2 - dy2 = zpl. 



Thus, our possible values for x appear as the first coordinate of one of 
the solutions of at most 2^'^ — 1 = 8191 Pell equations. For a given Pell 
equation, the sequence {yn/yi)n>i forms a Lucas sequence with real roots. 
The Primitive Divisor Theorem for Lucas sequences with real roots (see, for 
example, |3j, or the more general result from [1] which applies to all Lucas 
sequences) says that if n > 6, then y„ has a prime factor which is at least as 
large as n — 1. In particular, for the cases treated in [6j (and [7j), it suffices 
to check the first 42 (respectively 100) values of the component x of the Pell 
equations involved there and among these one finds all possible solutions of 
the equations considered there. 

We follow the same approach in the present paper. The first step is, 
as previously, to determine the form of all the possible solutions. This is 
done as explained above via the Primitive Divisor Theorem for the second 
component of the solutions to Pell equations. Since there are 25 primes p < 
100, this leads to first solving 2^^ — 1 Pell equations x'^ — dy^ = 1, the largest 
one of them having d with 37 decimal digits. This is clearly impossible 
using standard algorithms like the continued fractions, as the computations 
would be too slow and the fundamental solutions of some of the involved Pell 
equations would have hundreds of millions of decimal digits. Instead, the Pell 
equations are solved by first computing the regulator of the ring of integers 
of the corresponding quadratic field, and then from the regulator obtaining 
a compact representation of the fundamental solution of the Pell equation. 
The only algorithm known fast enough to compute the huge amount of 
regulators needed is Buchmann's subexponential algorithm. The output of 
this algorithm gives exactly the regulator under the Generalized Riemann 
Hypothesis and unconditionally is only a multiple of the regulator. Next, 
we check which of the solutions to the Pell equations will lead to solutions 
to equation ([1]). Finally, a check is performed that proves our search misses 
no solutions, thus eliminating the apparent dependence of our results on the 
Generalized Riemann Hypothesis. 
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2 Application of the Primitive Divisor Theorem 

Here, we explain in more detail the applicability of the Primitive Divisor 
Theorem alluded to in Section [1] to our problem. 



Let X be an integer such that x^ — 1 is a product of only the primes up 
to 97. We can then write x^ — 1 = dy'^, or, equivalently, x^ — dy'^ = 1, where 

d = 2"! • 3''2 . . . Qja25^ g. ^ |Q^ 1} for i = 1, . . . , 25. 

The only restriction for d above is that not all the ajS can be zero. This is 
a Pell equation, so (x,y) = (x„,y„), where Xn + y-avd = (xi + yivd)" for 
some positive integer n and xi + yiyd is the fundamental solution. As we 
have 2^^ — 1 possibilities for d, to get a finite number of possible solutions 
for our equation, we need an upper bound for n. 

Let d be fixed, rj = xi + yiVd and ( = xi — yiVd. It can be easily seen 
that 

'xn + ynVd = {xi + yiVdY = rf , 

^Xn-ynVd = (xi-yiVd)"- = C- 
From here, we deduce that 

Xn = :: and y„ 



2Vd 

It is easy to see that yi divides y„. We define m„. = Vn/yi- As Un = 
{rf- — C^)/{r] — C), it follows that (n„)„>o is a Lucas sequence of the first 
kind with real roots r] and C- By a result of Carmichael (see [3]) known 
as the Primitive Divisor Theorem, it follows that if n > 12, then Un has a 
primitive divisor p. This primitive divisor has several particular properties, 
the most important one for us being that it satisfies the congruence p = ±1 
(mod n). This implies that for n > 98, there exists a prime p > 101 dividing 
Un- Thus, n < 98. 

3 The algorithm 

We now give our algorithm. Let 

S = {2, 3, 5, 7, 11, 13, 17, 19, 23, 29, 31, 37, ... , 97}, 

be the set of all primes p < 100. The set S has 25 elements. Let Rd be 
the regulator of the ring of integers of the quadratic field Q(vw) and let 
xi + yiVd be the fundamental solution of the Pell equation x^ — dy"^ = 1. 

For all D G V{S) 
{ 

1- d = UpeDP 

2. Compute mRd 

3. Compute a compact representation of Xm + ymVd 

4. For i = 1 to 25 compute ordp-ym 
5_2 = 2°'''^pi^'" • . . . 97°''<ip25?^'" 



6. Compute all the convergents Pn/ln of the continued fraction expansion of 
\/d having q^ < z, and check whether p'^ — dq^ = 1 

7. If Rd - log 2 - log \/d Rs log z 

{ 

m = 1 and xi is a solution 

8. For i = 2 to 98 

{ 

For j = 1 to 25 compute oidp.yi 

z = 2°'''^pi?^» • . . . 97°'''ip25J'» 

\^ i ■ Rd — log 2 — log \/(f f« log 2, Xi is a solution 

} 

} 

} 

The algorithm searches through all 2^^ — 1 possible d. 

In step 1, a value for d is chosen. 

In step 2, Buchmann's subexponential algorithm is used to compute Rd- 
This algorithm returns a multiple mRd of the regulator Rd, unconditionally . 
If the Generalized Riemann Hypothesis is true, then m = 1. We will remove 
the dependence of our algorithm on the Generalized Riemann Hypothesis in 
step 6. 

A compact representation of an algebraic number /? £ Q(vd) is a repre- 
sentation of (3 of the form 

"^nli) . (3) 

where dj G Z, ctj = [aj + hj\/d)/2 G Q(v«), Oj, 6j G Z, j = 1, . . . , /c, and 
/c, a and dj have 0(log\/d) digits. 

In step 3, a compact representation oi Xm+VniVd is constructed using the 
algorithms described in [5]. The reason for using compact representations 
is that the standard representation of the fundamental solution of the Pell 
equation has 0(va) digits. Using the standard representation would make 
these computations impossible. More details about compact representations 
can be found in |5]. 

Since we only have the compact representation of the Xm + UrnVd, in 
step 4, we use algorithms from [9j to perform modular arithmetic on the 
compact representation. The p-adic valuations of ym when p is one of the 
first 25 primes are computed in the following way. We define v = Y[i=iPi- 
We first compute y™, (mod v). If y^ = (mod v), we conclude that ym is 
divisible by all 25 primes. Next, for all primes that we now know divide ym, 
we compute ym (mod pf), where c is a sufficiently large constant. We first 
take c = 15 and if y™ = (mod p^) still, we then replace c by 2c and repeat 
the computation. 



In step 5, z, which is defined to be the part of y^ divisible by the first 
25 primes, is computed. 

The purpose of step 6 is to remove the dependence of this algorithm 
on the Generalized Riemann Hypothesis. If the Generalized Riemann Hy- 
pothesis is false, it is possible that, without this check, we could miss some 
solutions in our search. Suppose therefore that m > 1 and that yk is such 
a solution that we are missing, meaning that y^ is divisible only by some of 
the first 25 primes for some k that is not a multiple of m. But yi | yt, so 
yi is divisible only by some of the first 25 primes. Also yi \ y^, meaning 
that yi divides the part of ym that is divisible by the first 25 primes, which 
in our notation is z. In other words, yi | z, implying yi < z. As yi has 
to be the denominator of a convergent of the continued fraction expansion 
of yd, it follows that by checking that the relation p^ — dq^ ^ 1 holds for 
all n satisfying g„ < z, we arrive at a contradiction. This proves that ei- 
ther m = 1, or that y^ has a prime factor larger than 100 for all positive 
integers k. This implies that our algorithm finds all solutions to equation 
([1]) unconditionally, and if ym = z, then m = 1. In practice, z will be a 
relatively small number, so usually only about 10 convergents will need to 
be computed. 

Since xi ss yi^/d, it follows that 

Rd = log(xi +yl^/d) ^ \og{2yiVd) = log 2 -^ log \/(i -h log yi , 

so we can determine whether yi = z by the test in step 7. In this test, we took 
that a Si 6 if I a — 6| < 0.5. If yi ^ z, then yi > 101 -z, so logyi > 4.61-1- log z. 
This shows that great numerical precision is not needed here. Just in case, 
we used 10 digits of precision. 

Step 8 checks whether any of the Xn are solutions for n = 2 to 98. With 
the purpose of speeding up the algorithm, the fundamental solution was not 
powered when computing ordp^yj. Instead, we computed the fundamental 
solution modulo the appropriate integer and then powering modulo that 
integer. Also, usually not all n have to be checked. This is a consequence 
of the fact that if a prime p divides yk, then it divides yik for every positive 
integer I. For example, if we get that 2/2 is divisible by a prime larger than 
100, then ^4, y%,... do not need to be checked. 

The running time of the algorithm is dominated by the computation of 
the regulator in step 2. Step 2 makes one the algorithm run in subexpo- 
nential time. It is the only part of the algorithm that is not polynomial. 
The computations were performed on a Intel Xeon E5430. The computation 
of the regulators took around 12 days of CPU time, while the rest of the 
computations took around 3 days of CPU time. 

Suppose that one wants to find the solutions of P{x^ — 1) < K using 
this algorithm. By the Prime Number Theorem, there are approximately 
K/logK primes up to X. This means that the algorithm will loop 2 ^^ 



times. The product of all these primes will be of size 0(e). This means that 
step 2 will run in 0{e^^^°^^). Thus, the time complexity of the algorithm 
is O {exp{K / log K + ^/K log K)). 

4 Results 

Theorem 1. a) The largest three solutions x of the equation P{x'^ — 1) < 

100 are 

'19182937474703818751, 

332110803172167361, 

99913980938200001. 

b) The largest solution x of P{x^ — 1) < 100 is x = 4217. 

c) The largest solution x of P{x^ — 1) < 100 is x = 68. 

d) The largest n such that P{x'^^ — 1) < 100 has an integer solution x > 1 
is n = 10, the solution being x = 2. 

d) The largest n such that P{x'^^ — 1) < 100 has an integer solution x > 2 
is n = 6, the solution being x = 6. 

e) The equation Pix^ — 1) < 100 has 16167 solutions. 

f) The greatest power n of the fundamental solution of the Pell equation 
{xi + yi\/d)" which leads to a solution of our problem is (2 + -v/3)-^^; 
i.e., n = 18 for d = 3. The case d = 3 also gives the most solutions, 
namely 10 of them. 

Proof: 

The proof is achieved via a computer search using the algorithm from Sec- 
tion [3l Part b) is proved by finding the largest square of all the x, c) by 
finding the largest a cube, etc. D 

The largest solution x has 20 decimal digits, the second largest has 18 
digits, followed by 5 solutions with 17 digits and 10 solutions with 16 digits. 
All of the mentioned large solutions are odd. The largest even solution x 
has 15 digits. 

Theorem 2. Write the equation ([1]) as 

x2-l = 2''i---97''25^ 

Then the following hold: 

a) The solution with the largest number ofai ^ is x = 9747977591754401. 
For this solution, 17 of the aiS are non-zero. 
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b) The solution with the largest X;?=i Oi is x = 19182937474703818751. 
For this solution, X]i=i ^i — ^7. 

c) The single largest Oj appearing among all solutions is ai = 27 and 
corresponds to the solution x = 4198129205249. 

Proof: 

Again, this is done via a computer search using the algorithm from Section 

El D 

Dabrowski [3] considered a similar problem, where the prime factors of 
x^ — 1 consist of the first k primes pi, . . . ,pk. He formulated the following 
conjecture. 

Conjecture 3. The Diophantine equation 

X -l=p^' ■■■p^" 
has exactly 28 solutions (x; ai, . . . , at) in positive integers, as follows: 

a) (3; 3), 

h) (5; 3,1), (7; 4,1), (17; 5, 2), 

c) (11; 3, 1, 1), (19; 3, 2, 1), (31; 6, 1, 1), (49; 5, 1, 2), (161; 6, 4, 1), 

d) (29; 3, 1, 1, 1), (41; 4, 1, 1, 1), (71; 4, 2, 1, 1), (251; 3, 2, 3, 1), (449; 7, 2, 2, 1), 
(4801; 7, 1,2, 4), (8749; 3, 7, 4,1), 

e) (769; 9, 1, 1, 1, 1), (881; 5, 2, 1, 2, 1), (1079; 4, 3, 1, 2, 1), (6049; 6, 3, 2, 1, 2), 
(19601; 5, 4, 2, 2, 2), 

/; (3431; 4, 1, 1, 3, 1, 1), (4159; 7, 3, 1, 1, 1, 1), (246401; 8, 6, 2, 1, 1,2), 

g) (1429; 3, 1, 1, 1, 1, 1, 1), (24751; 5, 2, 3, 1, 1, 1, 1), (388961; 6, 4, 1, 4, 1, 1, 1), 

h) (1267111; 4, 3, 1,1, 3, 1,1,2). 

The main result of [3] is that Conjecture [3] is true for fc < 5. From our 
data, we confirm Dabrowski's conjecture in a much wider range. 

Theorem 4. Conjecture\^ is true for k < 25. 

Proof: 

This is done by simply factoring all x^ — 1, where x is a solution of ([1]). D 

The next theorem follows also trivially from our results. Recall that a 
positive integer n is K-smooth if P{n) < K. In particular, the main result 
of our paper is the determination of all the 100-smooth positive integers of 
the form x^ — 1. 



Corollary 5. Let t he the largest odd solution, t = 19182937474703818751 
and s he the largest even solution, s = 473599589105798 of equation ([1]). 

a) The largest consecutive 100-smooth integers are x and x + 1 where 
X = {t — l)/2 is the largest solution of P{x{x + 1)) < 100. 

h) The largest consecutive even WO-sniooth integers are t — 1 and t + 1. 

c) The largest consecutive odd WO-smooth integers are s — 1 and s + 1. 

d) The largest triangular WO-smooth integer is {t^ — l)/8. 

As we mentioned in the Introduction, the problem of finding two consec- 
utive X-smooth integers was examined by Lehmer in |6j in the sixties. At 
that time, he was able to solve the above problem for the values K < 41. 
The advance of both computing power and theoretical arguments (namely, 
the compact represenations) allow us to solve the much harder problem of 
finding consecutive X-smooth integers for any K < 100. Note that, as was 
already mentioned at the end of Section [31 the difficulty of this problem 
grows exponentially with K. 

Our results can also be applied to finding k consecutive i^-smooth inte- 
gers, for any integer k. We obtain the following results. 

Corollary 6. The largest integer x satisfying 

P{x{x + l)...{x + n)) < 100, 
for a given n, are given in the following tahle: 



n 


X 


1 


9591468737351909375 


2 


407498958 


3 


97524 


4 


7565 


5 


7564 


6 


4896 


7 


4895 


8 


284 



Proof: 

To find k consecutive 100 smooth integers, we first create a list of all pairs of 
consecutive 100-smooth integers. Every odd solution of our starting prob- 
lem will give us one such pair. This is because if x is an odd solution, then 
(x — l)/2 and (x + l)/2 are consecutive 100-smooth integers. Once this list 
is created, we search for overlaps in these pairs and obtain our results. D 



In the recent paper [10], Shorey and Tijdeman proved several extensions 
of some irreducibility theorems due to Schur. The main results of [10] rely 
heavily on Lehmer's results from [6]. Thus, replacing, for example Lemma 
2.1 in |10j by our results, it is likely that the main results from ^Oj can be 
extended in a wider range. 

Remark. The tables produced by our computations can be found on 



the webpage |http://web. math.hr/~fnajrriari 
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